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Abstract 

We analyze the anomalies of superconducting state within a simple exactly 
solvable model of the pseudogap state, induced by fluctuations of "dielectric" 
short range order, for the model of the Fermi surface with "hot" patches. The 
analysis is performed for the arbitrary values of the correlation length ^corr of 
this short range order. It is shown that superconducting energy gap averaged 
over these fluctuations is non zero in a wide temperature range above Tc — the 
temperature of homogeneous superconducting transition. This follows from 
the absence of self averaging of the gap over the random field of fluctuations. 
For temperatures T > Tc superconductivity apparently appears in separate 
regions of space ("drops"). These effects become weaker for shorter correla- 
tion lengths (,corr and the region of "drops" on the phase diagram becomes 
narrower and disappears for ^corr 0, however, for the finite values of (,corr 
the complete self averaging is absent. 

PACS numbers: 74.20.Fg, 74.20.De 
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I. INTRODUCTION 



Among the anomalies of electronic properties of high - temperature superconducting 
cuprates (HTSC) especially interesting are the properties of the pseudogap state, observed 
in a wide region of the phase diagram 0,@]. From our point of view, the preferential scenario 
for the formation of the pseudogap state is based on the existence (mainly in the under- 
doped region of the phase diagram) of strong scattering of current carriers on well developed 
fluctuations of "dielectric" short range order (of antiferromagnetic (AFM) or charge density 
wave (CDW) type) 0. This scattering leads to essentially non - Fermi liquid like renormal- 
ization of electronic spectrum in certain parts of momentum space close to the Fermi surface 
and near the so called "hot" spots or "hot" (flat) patches 0. The preferential nature of 
"dielectric" and not of "superconducting" pseudogap formation follows from the number 
experiments, the appropriate discussion can be found in Ref. 0. 

The major part of existing theoretical papers is dealing with pseudogap formation and 
its influence on the properties of the system in normal (non superconducting) state, while 
only few are consider the anomalies of superconductivity in the pseudogap state In 
particular, in Ref. we have analyzed superconductivity in a simple exactly solvable model 
of the pseudogap state, based on the model Fermi surface (in two dimensions) with "hot" 
patches For the description of the pseudogap state we have used an exactly solvable 
model first developed for the one dimensional case in Ref. and for the limit of very large 
correlation lengths of "dielectric" short range order. It was shown that the superconducting 
energy gap averaged over the fluctuations of short range order is, in general, non zero also in 
the temperature region above the "mean-field" temperature of superconducting transition Tc, 
corresponding (according to Ref. 0) to the appearance of homogeneous superconducting 
state in a sample as a whole. This lead us to the conclusion that for temperatures 
T > Tc there appear superconducting "drops", which exist up to the temperature Tco of 
superconducting transition in the absence of the pseudogap of "dielectric" nature. This 
effect was attributed in Ref. to the absence of self averaging of superconducting order 
parameter (gap) in situation when the correlation length of short range order fluctuations 
is larger than the coherence length of superconductivity (the size of Cooper pairs). 

Under the assumption of self averaging superconducting energy gap the effects of finite- 
ness of the correlation length of short range order fluctuations were analyzed in Ref. where 
we have considered the pseudogap influence on Tc, calculated the temperature dependence 
of the gap for T < Tc, and derived the Ginzburg - Landau expansion for T ~ Tc. In this 
case we have used our "nearly exact" solution for the pseudogap state induced by Gaussian 
fluctuations of short range order, proposed first in Refs. for one dimensional case and 
generalized for two dimensions in Refs. ||T(]| , p]T| . Within this approach it seems rather difficult 
to find the solution without the assumption of the self averaging superconducting energy 
gap. Note that the problem of the presence or absence of self averaging of superconducting 
gap was not studied well enough in most of the previous work on disordered superconduc- 
tors. In most studies the self averaging property was just assumed on "physical grounds", 
with the reference on significantly different length scales on which the superconducting gap 
changes (coherence length ^o) and characteristic scales for the electronic subsystem (inter- 
atomic distance or the inverse Fermi momentum), e.g. in the impurity problem [p!2|-[I^, or 
the correlation length C,corr oi short range order in our problem p|,^,^), where it may be 
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expected that the complete self averaging of the superconducting energy gap appears for 
Ccorr ^0 lljH- We are unaware of any work, where the problem of self averageness was 
analyzed within some exactly solvable model of disorder. 

One of the main aims of the present work is to perform precisely this type of analysis 
within very simple (though probably not realistic enough) one dimensional model of the 
pseudogap state, induced by fluctuations of "dielectric" short range order with finite corre- 
lation length, proposed in a recent paper of Bartosch and Kopietz []T3[. An exact solution 



found in this work is very similar to that of our previous studies and allows to perform 
a complete analysis of the self averaging properties of superconducting energy gap for two 
dimensional model of "hot" patches, proposed in Refs. PJ^JTlH. Besides we shall study the 
full temperature dependences of superconducting energy gap for a superconductor with the 
pseudogap of "dielectric" nature. 



II. SIMPLIFIED MODEL OF THE PSEUDOGAP STATE. 



Let us consider an exactly solvable model of the pseudogap state proposed in Ref. |T5| us- 
ing slightly different approach. Consider an electron in one dimension in a periodic potential 
of the following form: 

V{x) = 2Dcos{Qx + (f)) (1) 

with Q = 2pp — k, where pp - is Fermi momentum, and k <^ pp - is some small deviation 
from the scattering vector Q Electronic spectrum close to the Fermi level is taken in 
the usual linearized form: 

^l=^p = Vp{\p\ -pp) ^p-2pp = -^p ("") 

6 = ^p-Q = -^p - vpk = -^p - r/ (2) 

where we introduced the notation rj = vpk {vp - Fermi velocity), which will be widely used 
in the following. The potential (|I]) can be rewritten as: 

V{x) = /)g«2p_FX-ifca; _|_ jj* ^-i2ppx+ikx 

where we have introduced the complex amplitude D De^'^ . 

The solution of this problem is elementary. In the "two-wave" approximation of the 
usual band theory the one-electron (normal) Green's function, corresponding to the diagonal 
transition p — >■ p, takes the following form (in Matsubara's representation): 

giiiisnpp) = r + -■ r^- r + - = 

^£^n - 6 ie + i + ri 



^This choice of the vector of AFM or CDW superstructure corresponds, in general, to the case of 
incommensurate ordering and fluctuations. 
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where in the last equahty we have introduced the following short notations: = ^, En = s. 
We can also introduce the non diagonal (anomalous) Green's function, corresponding to the 
Umklapp process p ^ p — Q: 

1 1 
gi2{i£nPP -Q) = -. -^D*- + ... = 

D* . , 

(5) 



Let now (|I|) be random field. Following Ref. |T^ we shall consider the very special model 
of disorder, when we shall assume the randomness of the scattering vector deviation k, with 
Lorentzian distribution 0: 

where k = and ^corr - is the correlation length of short range order. Phase in @) is 
also considered as random with uniform distribution between and 2tt: 

VM) = [^- O<0<2vr 
1 for other values 

Correlation function of random fields V{x) at different points can be calculated directly and 
is given by: 

< V{x)V{x') >= 2D^ cos[2pf{x - x')\ exp[-fi;|a; - (8) 

where the angular brackets denote averaging over (j^) and (^. The random field with 
precisely this form of pair correlator was treated first in Ref. ||16[, as well as in Refs. f^-Q, 



though in these works the Gaussian statistics of the random field was also assumed 0. In 
the model under consideration the random field Vix) is, in general, non Gaussian [|15|. The 
Fourier transformation of (P) takes the form of characteristic Lorentzian, which defines the 
effective interaction of an electron with fluctuations of short range order |]^: 

This type of effective interaction was introduced in all papers on the pseudogap of "dielectric" 
nature cited above. 

Green's functions averaged over the ensemble of random fields with distributions 
(^) and (1^) are calculated by elementary integration. The average value of the anomalous 
Green's function is simply zero (after the averaging over (^), which corresponds to the 



^In fact this corresponds to a specific model of phase fluctuations of the potential ([l|). 

^For the Gaussian random field all higher correlators of the random field V{x) are factorized "a'la 
Wick" into products of pair correlators (|8[). 
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absence of "dielectric" long range order. The average of the Green's function is obtained 
by term by term integration of the perturbation series @ with (|^), so that: 

+ ^ D* ^ D i D* ^ D ^ + ... = 

l^n - ip l^n + + f'VFi^ f'^n " ip l^n + ^p + 't'VpK lE^ - 



i£n + Cp + ivpn 



{iSn - Q{i£n + ^p + ivFn) -\D\ 



(10) 



This is precisely the exact solution for the Green's function given in Ref. |T5| . 

In the following we shall also consider the model with fluctuating amplitude D of the field 
(0), so that the appropriate Green's function can be obtained by the averaging of (|T^) with 
some amphtude distribution Vd{D). In particular, we can take the amplitude distribution 
in the Rayleigh form [|7|,|8|,P^ : 



2D ( 

V,(D) = (11) 

In this case the additional averaging of correlators (|]) and (H) leads to the simple replacement 
D — > W . The average Green's function of an electron in this case becomes: 



POO 

G{ienP) = / dDVoiD)- 
Jo 



i^n + Cp + ^"^F/t 



- ^p)(^^n + ^p + ivpK) - \D\ 
_^ iSn + ip + ivpf^ 
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^/-p-C ' '^P ' ""-^"^ (12) 

iien-Qiten + ^p + iVFf^)-CW^ ^ ' 

where W determines now the effective width of the pseudogap. In the limit of the large 
correlation length of fluctuations of (|Tp, i.e. for ^corr ^ cc (k — > 0), the solution (|T2|) 
coincides with that found earlier in Refs. ITU for the case of the Gaussian random field. For 



finite k this solution coincides with that proposed in Ref. |TT[ during the formal analysis of 
approximations, used in Refs. in the analysis of the general problem of an electron in the 
Gaussian random field with pair correlator given by (|]). In Refs. [[TTiplSl it was shown that 



the density of states, corresponding to the Green's function ([T^), possesses a characteristic 
smooth pseudogap in the vicinity of the Fermi level, and the values of the density of states 
are quantitatively very close [^,^,|r^ (practically for all energies in the incommensurate 
case) to the values obtained in Ref. P], as well as to the results of exact numerical calculation 
for the Gaussian random field, performed in Refs. [P"H|-p]|] 0. 

If the random field (|l|) is created by fluctuations of some kind of "dielectric" order param- 
eter (e.g. CDW or AFM), distribution ([TT]) may describe Gaussian fluctuations, existing at 



^Using the approach of Ref. |^ in the present model it is easy to find also the two particle Green's 
function and accordingly the frequency dependence of conductivity jl^. Unfortunately, the specific 
form of "disorder" (random field) leads to unphysical behavior at zero frequency, corresponding to 
an "ideal" conductor. 
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high enough temperatures |TO|)|ll| . For lower temperatures, even before the appropriate long 
range order appears, the amplitude fluctuations of the order parameter are "frozen out" (cf. 
PJ21[]) and we may assume the amplitude to be non random and put D = W, while phase 
fluctuations remain important up to much lower temperatures. Accordingly, the solution of 



the form of (|T0[), leading to sharp enough pseudogap for large correlation lengths ^corr |]T6 
can be used to describe the low temperature region of fluctuations of short range order. As 
we do not analyze the microscopies of "dielectric" fluctuations, all the parameters, charac- 
terizing these fluctuations, such as correlation length ^corr = and amplitudes D and W 
(the energy width of the pseudogap) are considered as phenomenological. The "low temper- 
ature" and "high temperature" regime of fluctuations of short range order can, in principle, 
be realized at different temperatures in comparison to the temperature of superconducting 
transition. 

Generalization to the case of two dimensional system of electrons, typical for HTSC - 
cuprates, can be done in the spirit of "hot patches" model for the Fermi surface, proposed 
in Refs. In this case we shall assume the existence of two independent types of fluctu- 

ations of the type of (|1|) [|, oriented along the orthogonal axes x and y, strongly interacting 
only with electrons from flat parts of the two dimensional Fermi surface, orthogonal to these 
axes. Accordingly, we assume the factorized form of two dimensional (random) potential 
for electrons: V{x,y) = V{x)V{y) l^j-Q. The size of these flat ("hot") parts of the Fermi 
surface is determined by an additional parameter a, so that 2a gives the angular size of the 
flat part, as seen from the center of the Brillouin zone 0]^-^. In particular, the value of 
a = ir/A corresponds to the square like Fermi surface (complete nesting), when all Fermi 
surface becomes "hot". For a < it/A the "cold" parts of the Fermi surface appear, where we 
neglect the scattering on fluctuations of the "dielectric" order parameter and electrons are 
considered as free. In this model all the physical characteristics, determined by the integrals 
over the Fermi surface, consist of additive contributions from "hot" and "cold" parts. Pseu- 
dogap renormalization of electronic spectrum takes place only on "hot" parts (and close to 
these parts), while on "cold" part the usual Fermi liquid (gas) behavior remains 0. 

This picture is in qualitative agreement with a numerous ARPES experiments on under- 
doped cuprates which show that pseudogap anomalies appear in the vicinity of (0, tt) 
point of the Brillouin zone, and vanish as we move to the diagonal direction. The presence of 
flat parts on the Fermi surface of HTSC - cuprates was also reliably observed in the ARPES 
experiments by several independent groups 0. 



III. GOR'KOV'S EQUATIONS IN THE PSEUDOGAP STATE. 

To study superconductivity for the system with pseudogap due to fluctuations of "di- 
electric" short range order we shall assume the simplest BCS form of pairing interaction, 
characterized by the attraction constant V, which is non zero in some energy interval 2uJc 



^Note the crude analogy of this picture to the concept of phase separation in HTSC - cuprates 
(stripes) ||2^, if we associate the correlation length S,corr with characteristic size (period) of stripes 
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in the vicinity of the Fermi level {uc - is the characteristic frequency of the quanta, re- 
sponsible for electron attraction). We have already used the same approach in Refs. 
In the present work we shall limit ourselves only to the analysis of s-wave pairing. There 
are no serious obstacles within our approach for the analysis of d-wave pairing, typical for 
HTSC - cuprates, though in this case the presence of the angular dependence (anisotropy) 
of superconducting gap leads @,|] to the additional integration over the angle and to a 
significant growth of time necessary for numerical computations. At the same time, it was 
shown in Refs. @-|^ , that the pseudogap influence on superconductivity is essentially similar 
both for s and ci-wave cases, differing only by the scale of parameters leading to the same 
changes of basic characteristics of superconducting state ((i-wave pairing is less stable to the 
dielectrization of electronic spectrum in comparison to s-wave pairing). 

On "cold" parts of the Fermi surface superconductivity is described by standard BCS 
equations. In the following we shall concentrate on the derivation of Gor'kov's equations 
for the one dimensional model, which is equivalent to the analysis of the situation at "hot" 
parts of the Fermi surface for two dimensional case [§,0. Green's functions (^, (|^) for one 
dimensional system in periodic field (|I]), can be written as matrix elements: 

iSn - 6 D* 
9ii = T- T\r- T-\ TT^ fi'12 



[ie^-ix){ie^-ix)-\D\^ {^en-il){^en-^2)-W ^ ' 

In the presence of Cooper pairing the Gor'kov's equations, constructed on "free" Green's 
functions of the type of (pTS]), are shown in diagrammatic form in Fig. |I|. In analytic form 
we have: 

Gil = gu - (711 AF+ - gi2^F^i 

G21 = 921 - g2i^F^i - g22^F^i 

F,\ = g;,A*Gn + gl^^^G^i (14) 

where superconducting energy gap is determined, as usual, from: 

/oo 
rfepi^iMp) = AT5:F+(£„) (15) 
~oo ^ 



np 



where we have introduced the dimensionless pairing constant A = Nq{{))V (Ai'o(O) - is the 
free electron density of states at the Fermi level.). 
The solution of Eqs. (|I4D gives: 



Gil = -;^[(^^ + ei)(g' + + D' + A') - D\i, + 6)] = 
^ {(z£ + Ok' + (e + ^)' + ^' + A2]+D'r/} 



Det 



F,\ = + ^2 + D' + A') = 

^ A*[e^ + {C + rif + D^ + A^] (16) 



Det 
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where 



Det = {e^ + ei +D^ + A^){e^ + + + A^) - + 6)'^' = 

= {e^ + e + D^ + ^^){e^ + {^ + vf + D^ + ^^)- V^D^ (17) 

where D denotes the real amphtude of fluctuation field In accordance with Eq. (|T3p 
Gor'kov's Green's function F^^ determines the superconducting energy gap. Taking into 
account the random nature of "dielectric" fluctuations Eq. ( [T3| ) should be averaged over 
the fluctuations of both "phase" r] = vpk and amplitude D, using distributions (H) and (for 
high temperature fluctuations) (0). 

Rather long, though direct, calculation of the integral in (|I5|) via residues gives: 
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ttA* 1 I e'^ + D^ + 



where we have introduced: 

e = Ve^ + A2 (19) 



Then from Eq. (|l^) we can immediately obtain the equation for superconducting energy 
gap in two dimensional "hot patches" model 



1 = 2ti\T I aJ^ie, A, r], D) + (20) 

n=0 ^ ^ 

where we have introduced the relative fraction of "hot" parts on the Fermi surface a = 
^a. The second term in ( PD| ) gives the standard BCS-like contribution from "cold" parts, 
occupying (1 — a) part of the Fermi surface. Summation over n in (^) is performed up to 
some maximal value determined by the integer part of the ratio 

Numerical solution of (^) allows to find the value of the gap A(?7, /)) for fixed 77 and D 
(i.e. for the given value of the random field of fluctuations (|l])) for any given temperature. 
After that we can perform averaging over (^) and (11) and find in this way temperature 
dependence of the average energy gap. In particular, for the "low temperature" regime 
of "dielectric" fluctuations it is sufficient to average only over the "phase" 77, so that the 
superconducting gap is given by: 



< A >= - r dv ,1^"", , A(r^, D) (21) 



In "high temperature" regime we have to add also the averaging over the amplitude D with 
distribution function (|TT|): 



< A >= rdDDe^p (-^] - r rfr^^^^A(r^,D) 



(22) 



As a result we shall find the temperature dependences of the average superconducting gap < 
A > without any statistical assumptions like the self averaging nature of the order parameter. 
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Analogously we can calculate the temperature dependence of dispersion < > — < A >2, 
allowing to estimate the randomness of A, i.e. the presence or absence of self averaging. 
Results of these calculations will be presented in the next section. 

As we have already noted in the Introduction, most papers on superconductivity in 
disordered systems actually assume the self averaging nature of the superconducting gap 
A. In this case A is treated, in fact, as non random and independent of parameters of the 
random field in which the electrons forming the Cooper pairs actually propagate. In our case 
these parameters are the amplitude D and "phase" t] of (0), accordingly the self averaging 
over these characteristics may be studied separately. 

Let A be self averaging over the fluctuations of rj. In this case we can treat A in Eq. (|I6|) 
is independent of rj. Then the anomalous Gor'kov's function averaged over the fluctuations 
of 7] can be written as: 

< >= ^ r dn^^ e^ + {^ + vY + D' + A^ , . 

71 J-oo ' 7]^ + v^K^e^ + + + A^)[e^ + + r]y + + A^] - r]W^ ^ ' 

This integral can be directly calculated, so that after long calculations we get: 

< F+ >= A* 



(l + ^) 




+e 






e 







(24) 



where e was introduced in (p^Q]). Accordingly we can calculate the integral of (|2^), entering 
the gap equation: 

roo 71 A* f 1 + ^ 

<FA>^ I d^< F+ >= , \ ' , (25) 



So, despite rather complicated form of the anomalous Green's function (123), in the gap 



equation the account of the interaction with fluctuations on "hot" (flat) parts of the Fermi 
surface reduces to more or less "standard" renormalization: 

e^e[l + -^] =£ 1+ ^ 



A^A(l + ^UAfl+ ] (26) 



analogous to that appearing in the problem of impurity scattering in superconductors p3 
Similar renormalization was already noted for the variant of the present problem in Ref . . 
The analogy with impurity problem here is rather natural as our parameter vpK = VFC,cJrr 
actually determines the characteristic inverse time of flight of an electron through the region 
of short range order of the size of ~ ^corr- Of course, the additional pseudogap influence is 
also connected with the appearance in Eqs. (p4|), ( p5|) of the square of dielectric gap D^. 

Finally, the gap equation determining superconductivity in the "hot patches" model with 
the assumption of self averaging over "phase" fluctuations takes the following form: 
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-1 _|_ VjpK -1 ~ 

1 = 2nXT { a ^ + 



2e 



+ £2 (l + 



(27) 



The solution of this equations, naturally, is simpler than that of Eq. ( pUD with later averaging 
over (pTI). In the absence of fluctuations of the amplitude of "dielectric" field D, which is 
valid for the "low temperature" region of fluctuations of short range order, it is Eq. (^) that 
determines "mean field" (in terms of Ref. |^) behavior of A(T) with respect to fluctuations 
of the random field (|l|). 

In "high temperature" region of fluctuations of short range order, assuming for D the dis- 
tribution function (0) and self averaging over the fluctuations of D, we obtain the following 
equation for the average superconducting gap: 

/n2\ i_l£fk 1 - n 

2-AT E ^ ^ /„ ^^^exp {-^] , +,2-- ^ + \ (28) 




describing the situation analogous to that studied in our previous work 0], where we have 
analyzed the influence of Gaussian fluctuations of "dielectric" short range order using the 
methods of Refs. |||,|^ . Here the fluctuations of (|l|) are treated exactly, but D is assumed to 
be self averaging. We shall see below that the results following from the solution of Eq. ( pH]) 
are very close to that obtained in Ref. 0. For n 0{^corr oo) Eq- (11) reduces to the 
similar "mean field" equation of Ref. [Q. The temperature of superconducting transition. 



determined by Eq. ( pT]) or Eq. (^), can apparently be identified as the temperature of 
the appearance of infinitesimally small superconducting gap, homogeneous over the whole 
sample 0. 

In the next section we shall present the results of numerical solution of Eqs. (P7|), ( P^D in 
comparison with the results of an exact analysis, based upon the approach using Eqs. (pOD, 

(0), (0)- 

IV. MAIN RESULTS AND DISCUSSION. 

Let us discuss the results of numerical analysis of gap equations derived in the previous 
section^. 

In Figs. D and ^ we show the dependences of the critical temperature Tc of supercon- 
ducting transition for the "low temperature" regime of "dielectric" fluctuations (at this 
temperature the "mean field" gap, determined by Eq. ( pTj) becomes zero) on the width 
of the pseudogap W (which in this case coincide with the amplitude of dielectric gap D) 
and correlation length of short range order. These results are in qualitative agreement with 
similar dependences for the "high temperature" regime of "dielectric" fluctuations (when 



^During our calculations we have assumed for the relative part of flat patches on the Fermi the 
value of d = 2/3 
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Tc is determined by Eq. (]28|)), as well as with dependences, obtained by us earlier for the 
different model of "dielectric" fluctuations with finite correlation length in Ref. §]. With the 
growth of the pseudogap width W the "mean field" value of Tc is suppressed. Diminishing 
correlation length leads to the "filling" of the pseudogap H,|8|,p!5|] and suppression becomes 
less effective. 

In Fig. ^ the curves show the temperature dependences of superconducting gap < A >, 
averaged over both amplitude D and "phase" rj ( "high temperature" region of fluctuations 



of short range order, where < A > is given by (p2|)), for different values of vpK. Dashed 
curves give the appropriate "mean field" temperature dependences of superconducting gap, 
obtained assuming the self averaging of superconducting order parameter over both the 
fluctuations of amplitude and phase, as determined by Eq. (p8|). 

Superconducting gap, averaged over the fluctuations of short range order, is non zero in 
the temperature region above T^, corresponding to the zero of the "mean field" gap (i.e. the 
gap homogeneous over the sample). More so, it is seen that the average gap is non zero also 
in a narrow region of temperatures larger than the transition temperature in the absence 
of short range order fluctuations (pseudogap) Tcq. This effect is due to the existence of 
fluctuations of the "phase" rj, when the Fermi level is in the energy interval, corresponding 
to the peaks of the density of states due to formation of dielectric gap. To understand this, 
note that for the given realization of "phase" rj and of amplitude D, the density of states 
has the form: 

1 ;,„5:,s(^^) JTPf^,^ '7 (29, 

iVo(O) rrNo{0) p [ for other values 

where gfi{Epp) - is the retarded Green's function, obtained from @) via standard analytical 
continuation ie„ + iO, and Nq{Q) - is the density of states at the Fermi level in 

the absence of short range order fluctuations. Thus for ^ D the Fermi level position 
approximately coincides with the peak of the density of states, leading to larger values of 
superconducting gap A{ri,D). The growth of the dielectric gap D leads also to the larger 
width of peaks in the density of states (p9D , so that for | ^ D superconducting gap A{r], D) 
grows with D. Then, for any temperature above Tco and for large enough amplitudes of 
dielectric gap D > D*[T), on the phase diagram in the variables r] and D there is always 
a narrow region close to the line | = D, where superconducting gap A[ri,D) is different 
from zero (cf. Fig.|^). This leads to the appearance in the temperature dependence of the 
averaged (over random field fluctuations) superconducting gap < A > of an infinite, though 
exponentially small "tail" in the temperature region above T^q. 

At the insert in Fig. ^ we show the temperature dependence of the relative mean square 
fluctuation of the superconducting gap SA/A = a/< A^ > — < A >^/ < A > for the "high 
temperature" regime of "dielectric" fluctuations. For large correlation lengths of short range 
order (^o/'Ccorr <^ 1) these fluctuations of superconducting order parameter are very strong 
for all temperatures, signifying the obvious absence of self averaging. Surprisingly, these fluc- 
tuations of superconducting energy gap are strong enough also for small enough correlation 
lengths, at least in the region of T > T^. In particular, the "tail" in the temperature depen- 
dence of < A > for T > Tc is observed even for vpn/Tco = 100, when ^o/C,corr ~ 30 ^ 1. 

The full curves in Fig. ^ show the temperature dependence of superconducting gap < 



A >, averaged over the "phase" rj (cf. (plD ), for the "low temperature regime" of "dielectric" 
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fluctuations, when the amphtude fluctuations of dielectric gap are frozen out and D = 
W. Dashed curves show the appropriate temperature dependences of the "mean fleld" 
superconducting gap, obtained with the assumption of self averaging superconducting order 
parameter over fluctuations of the "phase" r], and deflned by Eq. (^7]). For large enough 
correlation lengths of the short range order the average gap for T < Tc is very close to its 
"mean fleld" values and its "tail" in the region of T > Tc is relatively small. This behavior 
for the "low temperature" region of "dielectric" fluctuations is due to the fact, that for 
^corr there is no randomness in this model at all [i] = 0, D = W). Accordingly, 

the mean square fluctuation of the gap, shown at the insert in Fig. ^ is rather small for 
large correlation lengths and T < T^, but grows sharply for T > T^. As correlation length 
becomes smaller, fluctuations of superconducting gap 6 A for T < Tc at flrst grow , mainly 
due to the growth of randomness (parameter vpn determines the width of the distribution of 
fluctuations of 77), but afterwards diminish in the region of ^o/^corr ^ 1. In the "tail" region 
of the averaged superconducting gap (T > Tc) fluctuations of superconducting gap are quite 
large. Though these fluctuations diminish for smaller correlation lengths of the short range 
order C,corr, they still remain signiflcant even for small enough correlation lengths, i.e. for 

corr > 1. 

As well as in the "high temperature" region of "dielectric" fluctuations the "tail" in the 
temperature dependence of the average gap is observed here also in the region of T > T^o- 
It is explained by the same reasons as discussed above. However, for the "low temperature" 
region the amplitude of the dielectric gap is non random and flxed at D = W. Thus, for 
Tco < T < T*, where T* is deflned by D*(T*) = W, there exists a narrow region of "phases" 
close to ?7 = 2W, where the superconducting gap A{ri, W) is different from zero, while for 
T > T* such region is absent (cf. Fig.^. The value of T* determines the temperature up 
to which exists the "tail" of the average gap, i.e. the critical temperature for the average 
gap < A >. From the deflnition of T* it is obvious that this temperature is independent 
of correlation lengths and depends only on W. As the width of the peaks in the density of 
states (p9|), as well as A{ri,D), grows with the growth of D if condition ^ ^ D is fulflUed, 
the value of T* grows with W. The dependence of T* on W is shown at the insert in Fig. p. 



V. CONCLUSION 

In this work we have studied superconductivity within very simple model of the pseu- 
dogap in two dimensional electronic system which allows an exact solution. The central 
result is the explicit demonstration of the absence[21 e of complete self averaging of super- 
conducting order parameter (energy gap) over the random fleld of "dielectric" fluctuations, 
leading to the formation of the pseudogap. This is rather surprising from the point of view 



of the standard theory of superconductivity in disordered systems ||12Hl4|. The absence of 
self averaging is reflected in the appearance of fluctuations of the gap, especially strong for 
the temperatures larger than the "mean fleld" temperature of superconducting transition Tc, 
following from the standard equations, assuming the the self averaging nature of the order 
parameter. We identify this temperature with the critical temperature for the appearance 
of homogeneous superconducting state in the whole sample, while for the real disordered 
system the superconducting state is inhomogeneous and for T > Tc superconductivity exists 
in separate regions ("drops"), appearing due to random fluctuations of the local density of 



12 



electronic states. The difference with our previous work [Q, where this picture was discussed 
in the hmit of very large correlation lengths of short range order ^corr ~^ oo, the use of the 
model of Ref. [|I^ allowed to find the complete solution for arbitrary values of C,corr- This 
solution demonstrates the absence of complete self averaging of superconducting gap even 
for C,corr < C,o, which contradicts the naive expectations of the standard approach [Q. As 
we noted above, we are unaware of works, where the problem of self averaging of A was 
analyzed within any exactly solvable model of disorder. Here we presented precisely this 
type of analysis. It is unclear at present whether the results obtained will be conserved in 
more realistic models of disorder. 

It will be very interesting to analyze also the behavior of the spectral density and tun- 
neling density of states, analogous to that done in our previous work |^ in the limit of 
^corr ^ oo. It wiU bc especially important to consider the self averaging properties of the 
density of states, which is the standard assumption in the theory of disordered systems. 

As to comparison with experiments on high - temperature superconductors, we note Refs. 



2^ , |25[| , where the scanning tunneling microscopy (measuring the local density of states) on 
films of Bi2Sr2CaCu20s+5 clearly demonstrated the existence in this system of microscopic 
superconducting regions along with dominant semiconducting regions with typical pseudo- 
gap in the electronic spectrum. These observations are in qualitative agreement with the 
main conclusions of the present study. 

This work was partly supported by the grants 99-02-16285 of the Russian Foundation of 
Basic Research and CRDF No. REC-005, as well as by the program of fundamental research 
of the Presidium of the Russian Academy of Sciences "Quantum macrophysics" , and by the 
State Contracts of the Ministry of Industry and Science 108-11(00) (program "Statistical 
Physics") and 107-1(00) (program of HTSC research). 
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FIG. 1. Gor'kov's equations in one dimensional periodic field. 
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FIG. 2. Dependence of the critical temperature of superconducting transition for the "low 
temperature" region of "dielectric" fluctuations on the width of the pseudogap W for different 
values of correlation length of short range order =: 1. 0.1; 2. 1; 3. 10; 4. 100. 
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FIG. 3. Dependence of the critical temperature of superconducting transition for the "low 
temperature" region of "dielectric" fluctuations on the correlation length of short range order for 
different values of the pseudogap width ^ =: 1. 1; 2. 3; 3. 5; 4. 10. 



16 



2.0 



1.6 



0.8 



0.4 



0.0 









8 - 
^ 




ill 




4 




4 - 
2 - 




_^^^=============^=^^^^^-~~~~^^ 3__^ 

" ~ 








- 


r 


1 ' 1 ' 1 ' 1 ' 1 

0.2 0.4 0.6 0.8 l.O 
\ T/T 




3 










2 








^ \ 




1 


\ 


\ 

\ 










\ 

\ 

\ 

\ 

1 


\ 











0.0 0.2 0.4 0.6 0.8 1.0 1.2 



T/T 

FIG. 4. Temperature dependence of superconducting energy gap for the "high temperature" 
region of "dielectric" fluctuations. Full curves - superconducting gap < A > averaged over the 
amplitude D and over the "phase" r/, as determined by Eq. (22). Dashed curves - "mean field" 
superconducting gap determined by Eq. (28). At the insert - temperature dependence of the 
relative mean square fluctuation of superconducting gap. Curves are given for y~ = 3 and =: 
1. 0.1; 2. 1; 3. 10; 4. 100. 
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FIG. 5. The regions of the phase diagram with non zero superconducting gap for different 
temperatures above Tco, shown for the values of T/Tco- 1. 1.05; 2. 1.1; 3. 1.2. Dashed Hne 
corresponds to D = rj/2. 
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FIG. 6. Temperature dependence of superconducting gap for the "low temperature" region of 
"dielectric" fluctuations. Full curves - gaps A averaged over "phase" rj for fixed values of the 
amplitude D = W , as determined by Eq. (27). Dashed curves - "mean field" gap, determined by 
Eq. (27). At the insert A - temperature dependence of the relative fluctuation of superconducting 
gap. Curves are given for ^ = 3 and ^ =: 1. 0.1; 2. 1; 3. 10; 4. 100. At the insert B - 
dependence of the critical temperature T* on the width of the pseudogap. 
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